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By employing a mapping to classical anharmonic oscillators, we explore a class of solutions to the Nonlinear 
Schrodinger Equation (NLSE) in 1+1 dimensions and, by extension, asymptotically in general dimensions. We 
discuss a possible way for creating approximate rogue wave like solutions to the NLSE by truncating exact 
solutions at their nodes and stitching them with other solutions to the NLSE. The resulting waves are similar 
to compactons with the notable difference that they are not localized but rather extend over all of space. We 
discuss rogue waves in a 4 field theory in the context of a discretized Lagrangian and rogue wave behavior is 
shown to evolve into a steady state. Due to time-reversal invariance of this theory, the steady state found could 
alternatively evolve into a rogue wave giving rise to a large wave which seems to appear from nothing. 
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Introduction. The Nonlinear Schrodinger Equation (NLSE) 
is a very versatile equation used in many branches of physics, 
dictating the behavior of wave packets in weakly nonlinear 
media. It represents the evolution of optical waves in a nonlin- 
ear fiber OUT), the envelope of wave packets in ocean waves 
in an infinitely deep ocean t8lfT2ll. various biological systems, 
and the price of options in economics [13, 14 ]. One common 
way of solving the NLSE is through the Inverse Scattering 
Transform which uses the idea of Lax pairs. 0311161 

A very open problem associated with the NLSE is the de- 
velopment of rogue waves appearing in different media iTTvl— 
l20l . Rogue waves are large waves that seem to appear from 
nowhere and are, at least, 2-4 times larger in amplitude than 
surrounding waves. They are often preceded by a large de- 
pression referred to as a "hole in the sea" in oceanic terms. 
This is likely a nonlinear effect in most cases. 

Nonlinear Schrodinger Equation in 1+1 Dimensions. We 
consider the NLSE with a drift velocity and general nonlinear- 
ity in the presence of an external potential U and an additional 
source term g. In general dimensions, 



i—y = [-ic- 



dt 



dx 



-eV 2 ^- AF(|^| 2 )^ 
+U(r,t)9 + g(r,t)], 



(1) 



where e and A are tunable real parameters, F is a general func- 
tion of the squared modulus, and c is the velocity along the 
drift axis (chosen to be the x-axis). The general "source" term 
g does not make an appearance in usual NLSE. However, it 
can be treated following the transformation of the NLSE to 
a mechanical problem which underlies a central part of this 
work. We will ignore g and consider its effect only in a sepa- 
rate later discussion. In 1+1 dimensions, we choose the ansatz 



(2) 



with real u{x) > and e %b an arbitrary phase. This gives rise 
to an ordinary differential equation in u, 

d 2 u 

um = -e-—, - \F(u 2 )u + U(x, t)u. (3) 
dx 1 



Following a Galilean boost to the moving frame, x —> (x— ct) , 
and a subsequent interchange of space with time, Eq. [3] is the 
equation of motion for a nonlinear oscillator of spring stiffness 
uj subjected to an external linear force of strength U. Under 

,2 

this interchange, the term e corresponds to the inertial term 
of the mass times the acceleration, and [uju + \F(u 2 )u] is an 
effective internal forcing term. In this form, the strength of 
the nonlinearity is proportional to A. In what follows, we first 
consider the case of U = and then comment on non-zero 
U. Equation [3] is that of a classical particle in an effective 
potential given by 



Veff(u) = -U 2 + \ 



/ dv[F(v 2 )v]. 



(4) 



The energy of a classical oscillator is a constant of motion, 
this case, the corresponding "energy" is given by 



In 



E = V eff (u) 



e du 2 
2 W 



(5) 



A particular case is that of solitons for which u — » at large 
\x — ct\ (or large times in the corresponding mechanical prob- 
lem) and consequently E = 0. General periodic solutions 
appear for general non-zero E. 
Rewriting Equation [3] we have 

dV =-(6) 



dx = ± 



2E 



A d(v' 2 ) 



dx 



F(v' 2 )dv' 



For definitiveness, we consider the cubic nonlinearity 
F(u 2 ) = u 2 in our Eq.([6]). The effective potential is then 
V e ff = uju 2 /2 + A^ 4 /4. This effective potential is plotted in 
Figure [T] for uo = ±2 and A = 1. In general, when uj < 0, 
we have a double minimum of V e ff(u). In such cases, when 
E < 0, u is restricted to be positive definite. When uj > 0, we 
have a single miniumum at the origin (u = 0). In all cases, as 
we increase the value of E, more values of u become possible 
(E > V e ff(u)). A transition in the form of the effective po- 
tential (from a double well form to a single minimum) occurs 
at uo = 0. In passing, we note that this is similar to what oc- 
curs in Ginzburg-Landau theories when the coefficient of the 
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FIG. 1 . Plot of the effective potential shown in Equation[4]for a cubic 
non-linearity (F — u 2 ) for uj — 2 (single minimum) and uj — — 2 
(double minimum) with A = 1. 

quadratic term (in this case uj) switches sign at the transition. 
We find that in the case of F(u 2 ) = u 2 , 

^ V4£A + uj 2 

" = -y-A + — a — x 

cn J ^AT^ 1_ uj 

where cn is a Jacobi Elliptic Function. This is consistent with 
companion works l2Tll22ll . The Jacobi Elliptic Functions are 
doubly periodic in (x — ct). In this case, the periodicity of 
the Jacobi Elliptic integral allows a harmonic (Fourier series) 
expansion of u l23ll . 

NLSE in d+1 dimensions. We consider the NLSE of Eq.Q 
in d+1 space-time dimensions. With r = \r\, r = rn, and 
angular momentum operator, L, in d spatial dimensions, the 
Laplacian becomes V 2 = + ^£ - ^. If L 2 ^ is 
bounded in the co-moving frame at speed c, for asymptotically 
large r along any direction n on the unit sphere, Eq.Q reads 

9 T / * n r ^ 2 t / * \ 
i — ^{r,n) = [- e _*(r,n) 

-AF(|^| 2 )^(r,n) + /7(r,t)^(r,n)]. (8) 

Along any direction n, self-consistent solutions are enabled 
by the solution of the 1+1 dimensional NLSE with the cor- 
responding potential U. Here,"self-consistency" is meant to 
imply the resulting wavefunction ^ has dependence on n such 
that L 2 \I/(r, h) is bounded. The irrelevance, at large r, of the 
two terms ^— ^ §^ and for bounded finite wavefunctions 

r or r z 

^ implies all large r solutions are those of the I = (s-wave) 
type with a radial dependence identical to that 1+1 dimen- 
sional systems investigated in Eq.([6]). A special subclass of 
bounded wavefunctions of the form of Eqs.Q [6]) are those 
which vanish as r — >• oo, of pertinence to many usual (linear 
Schrodinger type) problems for which j d d r\^f\ 2 is normal- 
ized. Such a vanishing of the wavefunctions renders the higher 
order non-linearity irrelevant at large r and reduces the prob- 
lem to the usual linear Schrodinger equation. We consider 
more general waves (idealized deep ocean type and others) 
for which the wavefunctions need not vanish at spatial infin- 
ity. In this simplest s-wave type case, the behavior along all 
directions n is the same. In 2+1 dimensions, with an ansatz of 
the form of Eq.j^, i.e., V(x,i) = e- iujt (Y,m u m(r)e irncf> )e i6 



with the cylindrical coordinates (r, <fi) taken in the moving 
frame at velocity c along the drift (x-) axis, the NLSE in the 
case of cubic non-linearity reads 

/ •> r/ d 2 Id. 
H~A ^ ^ u rni u rn2 u rn — rni -^ rn2 ]. (9) 

mi ,rri2 

In the large r limit, under the exchange of the radial coor- 
dinate with time, Eq.([9]), for all bounded wavefunctions with 
finite radial derivatives, is that of classical oscillators with dis- 
placements {u m } with a uniform non-linear coupling A. 

Noise. For a non-zero U in Eq.Q that depends only on the 
coordinates in the drifting frame at speed c, the mechanical 
analogue problem of Eq.([3]) along the radial direction (or, in 
general, in 1+1 or 2+1 dimensions) is that of an anharmonic 
oscillator subjected to a linear force uU (t). Anderson local- 
ization [24-28] will occur for random U (t). This means that 
an instanton-like solution of the mechnical problem is non- 
zero only in a finite interval of time. This corresponds to a 
spatial localization of the wavefunction ^ over the finite inter- 
val of \r — cte x | with e x a unit vector along the drift direction. 

Next, we briefly consider the case of a periodic potential U 
of period a for A = 0. The invariance of the second order 
differential equation of motion of Eq.Q under translations by 
a implies that independent linear solutions scale as /n x ^ a U(x) 
with general complex \i which may be of unit modulus or real. 
The case with complex |/i| = 1 constitutes a mechanical ana- 
logue of Bloch's theorem. Normalization demands that for an 
ideal periodic potential, |/i| = 1. For real /i with > 1, 
there is an instability (parametric resonance) wherein u ex- 
hibits unbounded increase. Such a case indeed occurs for 
a linear oscillator when the potential is U = h cos jx with 
1 7 — 2^/uj/e\ < h^/uj/Ae wherein when coupled to the unper- 
turbed u, the noise term uU acts as a periodic external force 
with a Fourier component having a period that is close to that 
of the natural resonant period of the unperturbed oscillator. 
This leads to large oscillations. A similar resonant Fourier 
component is generated by a cubic term in u. It is conceivable 
that such a situation is emulated for nearly periodic U over 
some spatial range also in the presence of non-linearities. We 
speculate that noise and effects from the ocean floor along 
with collisions with other waves may act as effective bound- 
ary conditions on the waves similar to those discussed in l29l 
which may trigger wave amplification in some cases. 

Source terms. Following the outlined prescription of con- 
verting the NLSE with the ansatz of Eq.Q into a classical me- 
chanical problem where the spatial coordinate of the NLSE 
in the moving frame is replaced by a time coordinate, / be- 
comes an external force. In particular, if in the 1+1 dimen- 
sional rendition of Eq.|l]) we have g = e luJt f(x — ct), then 
instead of Eq.([3]) we will have, in the moving frame at speed c, 
uju = —e^i — \F(u 2 )u + Uu + f(x). Under the interchange 
of space with time, this corresponds to a classical non-linear 
oscillator subjected to external forces (Uu + /). The results 
of the classical analysis are then replicated. For instance, in 
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the case of U = A = and an oscillatory external "force" 
f = fo cos fit, a resonance appears when Q = \Jujje wherein 
u ~ t sm(flt + (j>). The solution of the homogeneous equation 
(with a vanishing U = f = 0) is that of Eq.([6]). Whenever the 
"force" / is of the form of the homogeneous solution, a reso- 
nance appears. In the case of a cubic non-linearity (A / 0), a 
resonance occurs when / is of the form of Eq.Q. 

Rogue Wave Type Solutions. We consider emulating rogue 
waves by spatially stitching solutions together at common 
nodes. Following the formalism associated with compactons 
1301 , we consider wavefunctions which are piecewise continu- 
ous with cutoffs at the nodes. We can truncate the function, u, 
at the nodes and append on another function of our choosing 
with minimal error as long as the value of the function and its 
derivative at these points is close to zero. This process differs 
from the process used to create compactons in that the wave- 
function extends over all space and is not confined to a local- 
ized region. The process that we briefly describe produces an 
approximate solution to the NLSE (including, as a particular 
case, the linear Schrodinger equation) with deviations from a 
solution of the form of Eq.Q only in the vicinity of the nodes. 
The initial form of these solutions can be expressed as 



#(a?,0) 




(10) 



where the common nodes of the function U2 and u\ and us 
are, respectively, denoted by x\ and X2- In higher dimen- 
sions, {xi } denote the coordinates of planes (either Euclidean, 
spherical or other) along which these solutions are stitched. 
Here, {u^} correspond to our solutions of the form of Eqs.|6j 
[7]). We note that perfect nodes where the function and its 
derivative are exactly zero cannot exist, because for u of the 
form given in Eq.Q, the derivative is non-zero at all zeros of 
the function. Below, we once again set e = A = 1 and con- 
sider the cubic NLSE in the absence of an external potential. 




FIG. 2. Plots of the function u at t — for A 
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As uj becomes increasingly negative, the oscillations be- 
come narrower and taller as in Figure [2(a)! For negative u, 
negative E is allowed so long as E > f. A form for neg- 



ative E is shown in Figure 2(b) Although wavefunctions \I/ 
formed by piecewise stitching such solutions are special, it is 
clear that in order to have transient rogue wave type phenom- 
ena, the current j cannot be spatially uniform (with a velocity 




FIG. 3. Plot of |^| of Eq. [I0|with c = 1, u 2 = -10, ui = uj 3 = 
-1, E = lO" 12 , and x\ = 5.5 and x 2 = 20. In this plot, we shift 

(x — ct — xo) — >> x with xq — 15. 



c) in the rogue wave region. Far away where the current is of 
uniform speed c, it is of the form given by Eqs.([6|[7]). 

A wave U2 of large amplitude by comparison to wi and 
^3 emulates a rogue wave. A plot of ^(x : t — 0) is pro- 
vided in Figure pi The continuity equation |^ + V • j = 
with the density given by p = |\I/| 2 and the current j = 
ie(^W* - \E*V\EQ applies to the NLSE. For wavefunctions 
of the form of Eq.( 10) with real uu the current j = at all 
points apart from the nodes. If an interpolating form is chosen 
near the nodes, the current j assumes non-zero values only in 
the regions close to the nodes. If we consider the co-moving 
frame and define P = f* 2 pdx where p = |^| 2 , then integrat- 
ing the continuity equation gives 



dP 

~dt 



J1-J2 



(ID 



where ji andj2 are the currents at x\ and X2, respectively. If 
we consider a rogue wave centered about x = extending 
from x = —a to x = a and define P Q = j^^pdx with 
X ^> 0, then with P the integral over the region of extension 
of the rogue wave, we can define the lifetime r of the rogue 
wave as the time for P in the region of the wave to arrive at 
P n . We then find that 



dP 



p 32 ~ Ji 



(12) 



Alternatively, to avoid the possible problem of slow conver- 
gence to P Q , we can consider the half life of the rogue wave: 



TL/2 



dP 

(P+P )/2 32 - jl 



(13) 



We can concretely see the behavior of a rogue wave type 
solution by evolving rogue wave initial conditions according 
to the linear, free-particle Schrodinger equation (LSE). The 
behavior of wavefunctions under the LSE should give us qual- 
itative understanding of the time evolution of rogue waves. 

In the co-moving frame, we can write down rather gen- 
eral initial conditions for a rogue wave as ip(x, 0) = cos x + 
ae~ x , with (1 + a) 2 the amplitude of the rogue wave rela- 
tive to the background, typically between 2 and 4. Evolving 



4 



i/;(x, 0) according to the LSE, we find that 



2 cos x + a : 
VI + lit 

The amplitude of the rogue wave at the origin is 

a 2 vT 



2it 



Vl + 4t 2 Vl + 4t 2 



(14) 



c.c. (15) 



After an initial drop in amplitude, the amplitude of the oscil- 
lations at x = decays to zero as ^ . 

77ze AfL&E and Rogue Waves in a (j) 4 Model. We next con- 
sider the (j) 4 theory given by the Lagrangian density 

.2 



c = 



1 



A 



(16) 



m 

2 - . ~Y^- 4V 

with a complex scalar field and implicit summation over 
space-time coordinates fi. The Euler-Lagrange equation, in 
the non-relativistic limit, is a NLSE with a potential shift, m: 

= V 2 | A 
2m 12m 

Solutions are given by those of Eq. |l] ) with F = | ^ | 2 mul- 
tiplied by e~ zmt . If we discretize the <p 4 Lagrangian density 
and consider a 2-dimensional lattice of points with spacing A, 
the imaginary time (t = it) Lagrangian becomes 

m 2 A 2 



mo 



(17) 



i,3 



'■Mi. 
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2A 2 



13 

A 2 ^ 



(18) 



where the indices i, j denote the two Cartesian coordinates of 
points in the plane, a denotes each pair (i, j), (fc, Z), and £ a is 
zero for non-adjacent and one for adjacent oscillator pairs. 

In the linear (A = 0) case, the system is quadratic and solv- 
able by diagonalization. We analyzed the time evolution for 
different parameters by a perturbative scheme in A. In order 
to underscore the ubiquitous nature of rogue waves (includ- 



ing, as noted in some of our earlier discussions (e.g., Eq.( 10), 
their presence already at the linear level), we plot below the 
evolution of an initial rogue wave state for the coupled linear 



oscillator Lagrangian of Eq. ([18]) with the parameters m = 2 
and A = 1. The even appearance and low amplitude of the 
surface at large times illustrates that, already at a linear level 
(i.e., that with A = 0), the dispersion of the normal modes of 
the unperturbed sysetm (A = 0) enables an evolution of steady 
states into rogue wave configurations. As the theory is time- 
reversal invariant, we can conclude that a steady state given 



by Fig 4(b) could evolve into the rogue wave in Fig. 4(a) 



Conclusions. We investigated the NLSE by considering a 
corresponding classical mechanics problem of non-linear os- 
cillators. The focus was on determining oscillatory waves 
with a general drift velocity c. By invoking this analogy, 
(1) We determined solutions to general NLSE with a gen- 
eral drift velocity c in 1+1 dimensions for an arbitrary non- 
linearity F in Eq.Q . The general result is that of Eq.([6]). 



\<M 2 2 \ 




(a) 



(b) 



FIG. 4. (a) shows a region of the cell surrounding an initial Gaussian 
rogue wave on a cosine background where (i, j) denote Cartesian co- 
ordinates and 1 0| 2 denotes the modulus of the squared amplitude, (b) 
shows the same system at large times (t = 60, 000 iteration steps). 



(2) We discussed solutions in general dimensions for large ra- 
dial coordinates. 

(3) We discussed the effects of noise and external sources. 

(4) We demonstrated that solutions could potentially be trun- 
cated and attached to each other (as in Eq. ( fT0| ) in order to cre- 
ate approximate solutions that behave very much like a rogue 
wave. We thus introduced the notion of "generalized com- 
pactons " which unlike usual compactons (and solitons) have 
a non-zero (periodic) background. The lifetime of such solu- 
tions scales as the reciprocal of difference in currents in the 



vicinity of the stitching points as given by Eq.( 12). We qual- 
itatively saw the behavior of the amplitude of a rogue wave 
over time by solving the linear Schrodinger equation for rogue 
wave-like initial conditions. 

(5) We, specifically, employed considerations of the continu- 
ity equation to discuss life-times of rogue wave type states. 

(6) We reported on the appearance of rogue waves in <p 4 the- 
ories. Specifically, we modeled rogue waves in a system of 
coupled oscillators. We showed that this approach allows ini- 
tial rogue wave-like conditions to evolve into what appears to 
be a steady state. By time-reversal invariance of the equations 
of motion, we can conclude that rogue waves can evolve out 
of a steady state of that form. 
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